Weak localization in multiterminal networks of diffusive wires 
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We study the quantum transport through networks of diffusive wires connected to reservoirs 
in the Landauer-Biittiker formalism. The elements of the conductance matrix are computed by 
the diagrammatic method. We recover the combination of classical resistances and obtain the 
weak localization corrections. For arbitrary networks, we show how the cooperon must be properly 
weighted over the different wires. Its nonlocality is clearly analyzed. We predict a new geometrical 
effect that may change the sign of the weak localization correction in multiterminal geometries. 
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How to increase transmission with weak localization ? 
This is one the question addressed in this letter devoted 
to a general description of transport in networks of quasi- 
one-dimensional weakly disordered wires. At the classical 
level, the network is equivalent to a network of classical 
resistances, which gives the dominant contribution to the 
conductances. Additionally, there exists a small correc- 
tion, the "weak localization" correction, originating from 
quantum interferences. For the conductivity, it reads 
[1,2] : Acr(f) = -^Pc{f,r) (with h = 1, without spin 
factor), where the cooperon P^, is the contribution from 
pairs of time reversed trajectories to the return prob- 
ability. Networks are particularly well suited to study 
interference effects and the first experimental studies of 
weak localization in such systems, performed on honey- 
comb lattices [3], showed the oscillations predicted by 
Altshuler-Aronov-Spivak (AAS) [4]. These experimental 
results were well fitted by the theory of Dougot & Ram- 
mal (DR) [5] whose starting point is a uniform integra- 
tion of Aa{r) over the wires of the network [6]. However, 
although this procedure is meaningful for a translation 
invariant system, it is not valid in general for networks. 
In this letter we demonstrate that the correct expression 
of the weak localization correction to the resistance TZ of 



that a uniform integration of the cooperon is applicable 
only to regular networks in which the weights of the wires 
dTZ'^^ /dR'^^^ are all equal, like in the experiments of Ref. 
[3]. We show that in a multiterminal geometry, some 
of these weights can change in sign. This can lead to a 
change in sign of the weak localization correction. 




a network of wires {iiiy) of lengths 
instead : 



FIG. 1. A network of diffusive wires. The network is con- 
nected at the vertices a', 13' and ji' to external reservoirs 
through which some current is injected in the network. 

As on figure 1, we consider networks that can be con- 
nected to external contacts (here the contacts a', (3' 
and /i'). The transmission probability between two con- 
tacts (the conductance matrix element, up to a factor 
e^//i), when averaged over the disorder, can be written 

• •, where the first term is 



and section s is {Ta'fi') — T^mi + AT^/^ 



—^AR^, with 

(All/) 



= Idf- 



-Acr(f) 



(1) 



where TZ'^^ is the classical resistance obtained from the 
classical laws of combination of resistances i?^',^. The 
sum runs over all wires (/it^). (Tq is the Drude conduc- 
tivity. This result has a simple structure : it could be 
obtained from small variations AR^i, in the classical ex- 
pression TZ'^^ of the resistance. However it is highly non 
trivial since, due to nonlocality, it is not possible to get 
a quantum formula for the resistance of the network as a 
function of quantum resistances of the wires. It is in fact 
not even possible to define a quantum resistance for a 
wire, independently of the whole network. Eq. (1) shows 



the classical result (Drude conductance) and the second is 
the weak localization correction. Our aim is to give a sys- 
tematic way to compute the weak localization contribu- 
tion in terms of matrices encoding the information on the 
network (topology, lengths of the wires, magnetic fluxes). 
The most natural approach to describe transport in net- 
works is the Landauer-Biittiker approach. However some 
difhculties related to the question of current conservation 
are more conveniently overcome in the Kubo formula- 
tion. We first recall some features of the transport theory 
of weakly disordered metals in the Kubo approach and 
eventually use the connection to the Landauer-Biittiker 
formalism that we apply to the case of networks. Finally 
we consider several examples. 

Classical transport. The classical transport is described 
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by two contributions : a Dnidc contribution (short range) 
and the contribution from the diffuson (ladder diagrams) 
which is long range. Then, in the diffusion approxima- 
tion, we have [9] : 




((Ty (f,f '))class 



= <7o [5ij S{f- f') - ViV'jPdif, f')] = (70 cj^ij (r, f ') , (2) 

where the diffuson Pd is solution of the equation 

—APd{f,r') = S{f — r'). An important requirement 
of a transport theory is to satisfy current conservation 
Vicrij(f, r ') = 0, what the classical conductivity (2) does. 
The question of boundary conditions is not a trivial one. 
In most works, it is argued, hj analogy with the problem 
of classical diffusion, that n- VPd|ax> = on the reflecting 
boundaries of the domain "D (n is the vector normal to 
the surface), while Pd\dTi = on the boundaries where 
the disordered system is connected to a reservoir. This 
latter is described by a region free of disorder. However 
it was shown in [10-12] that the diffuson vanishes at a 
distance of order £e inside the reservoir, is the elastic 
mean free path. In the case of a quasi-ld wire, the dif- 
fuson vanishes at a distance Xd = where d is the 
dimension with ai = 2, = 7r/2, = 4/3. 
Current conserving weak localization correction. The 
weak localization correction involves the cooperon which, 
in a magnetic field B = V x A, is solution of : 



V - 2\eA 



Pc{f,f') = 6{r-r'). (3) 



The contribution of the cooperon reads [1,2] : 



(c7y(f,r')) 



cooperon 




^%5(r-f')Pc(r,rl. 



(4) 



It is clear that the additional contribution of the cooperon 
(4) does not respect current conservation. In the same 
way that the classical conductivity is built of short range 
(Drude) and long range (diffuson) contributions, the 
weak localization correction contains long range terms 
additionally to the short range contribution (4). A pos- 
sible way to build a long range diagram is to dress the 
current vertices with diffusons, like on figure 2. 




The connection of the cooperon and the diffusons is re- 
alized by introducing a Hikami box [13]. Dressing either 
the left current line with a diffuson, or the right one or 
both, we obtain three diagrams that add to the contribu- 
tion of the cooperon. However these dominant diagrams 
are not the only ones needed to satisfy current conser- 
vation [11]. This is related to the fact that the diagram 
2, which seems at first sight to be the only long range 
diagram needed to compute the weak localization, gives 
a volumic divergence V^(0) to the conductance (for de- 
tails, see [14]). This problem was first mentioned in the 
study of conductivity fluctuations [9] where a more sim- 
ple procedure to construct a current conserving theory 
was proposed, that avoids to consider all the set of dia- 
grams of [11]. Kane, Serota and Lee showed, for the con- 
ductivity fluctuations, that the expression of the current 
conserving conductivity is obtained by a "convolution" 
of the short range object with the function involved 
in (2). Finally the expression of the weak localization 
correction to the nonlocal conductivity reads : 

{Aaij{f,r')) = J dp dp' 

<l}ii'{r, p) 4>jj' (r", p') {ai'.j'{p, p')) 

cooperon ; (5) 

which obviously respects current conservation. 

Networks. We now consider specifically the case of net- 
works such as the one on figure 1. The transmission Ta'fj' 
between the two contacts (conductance in unit /h) is 
related to the nonlocal conductivity (Jij{f, r') with f be- 
ing integrated through the section of the contact a' and 
r ' through the section of the contact /3' [15]. For quasi-ld 
wires, we obtain the expressions 



Otd^c 



(6) 



where Nc is the number of channels in the wires, and 
2 /• 



I dx^Pd{a:,x) Pc{x,x) ^Pa{x,§;) , 

J Network '^^ 

(7) 



FIG. 2. A long range diagram for the weak localization. 



that involve the one-dimensional diffuson and cooperon, 
solutions of the diffusion equation [7 — D|] P{x,x') = 
S{x — x'), where D^; = d^^ — 2ieA{x) is the covariant 
derivative and 7 = l/L"^ (for Pd we set 7 = and 
A{x) = 0). The notation Pdi^',x) means that the diffu- 
son is taken at a distance ie of the vertex a'. Eq. (6,7) 
assume Dirichlet boundary conditions at the vertices con- 
nected to reservoirs. Indeed, it can be shown that the 
correct boundary conditions mentioned above can be re- 
placed by Pd\d-D = 0, providing to insert a multiplicative 
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factor 4, independent of d, as it has been done in (6,7). 
This simphfies shghtly the calculations. 

As a check, let us consider the case of a quasi-ld 
wire of length L connected to reservoirs at both sides. 
We recover the Drude dimensionless conductance — 
oidNcie/L, and the correction /S.g = — ( coth — 
which interpolates between At; = —1/3 in the fully co- 
herent limit [L <^ L^), and Ag = —L^/L for L Lip. 

We introduce the adjacency matrix UafS that encodes 
the information on the topology of the network : aap = 1 
if the vertices a and f3 are connected by a wire (a/3), 
whose length is denoted by lajs, cLap = otherwise. To 
construct the solution of the diffusion equation, we need 
to specify boundary conditions at the vertices. We im- 
pose continuity of P and '^a/3^(a/3)(Q^) = ^aP(pi) ■ 
the adjacency matrix constrains the sum to run over the 
neighbouring vertices of a. The parameters describe 
how the network is connected. Aq = for an internal 
vertex and Aq = cxd at the vertices connected to external 
reservoirs, which imposes a Dirichlet boundary condition. 
The solution for Pc involves the matrix [5,7,8,14] : 

Maf3 = Sa(3 ^Xa + y/^^Uaf, COth(^;Q^)^ 



-a.af3 



(8) 



The diffuson is expressed in terms of the same matrix 
with 7 = and no magnetic flux 6af3 = : 



1 



(9) 



This matrix encodes the information about the classical 
conductances aaNcie/l^iv of each wire {^v). Note that 



Aq' = oo (at a reservoir) implies that {M~ 
fi. The classical conductance is given by 



T3p> - ?^ (-Mo-^^ 



Q/3 



0, V 



(10) 



This result is only valid for a' ^ (3' [11]. It coincides with 
the one obtained for a network of classical resistances, as 
it should. 

In Eq. (7), we separate the integral over the network 
into contributions of the different wires : S^f,i-yf,aT:\i ~* 
I{fj.u)- Si^ce the diffuson Pdia',x) is linear in x 
on the wire (/iJ^), the derivatives produce coefficients de- 
pending on (/ii/) : 

X ~ {Mo^) + {Mo') + S^^o^S^c.' Lc.'] ■ (11) 

Replacing (11) in (7) shows explicitly the non trivial 
weights that should be attributed to each wire when inte- 
grating the cooperon over the wires. These weights have 
a clear physical meaning : they can be related to deriva- 
tives of the corresponding classical conductance with re- 
spect to the lengths of the wires. From (7,10,11) we ob- 
tain : 



a'l3' 



adNdc 



E 



dlau 



AxPc{x,x). (12) 



We have demonstrated (1), since ctq — ^-^^j^.The in- 
tegral of the cooperon over the bond {^v) is a nonlocal 
quantity that carries information on the whole structure 
of the network through the matrix Ai~' : 



AxPc{x,x) = TT-^l 



iM-') +{M-') 1 (cothV^;^.--^^^ 



sh ^l^,u 



-1 + coihy/^l^y ' 

sh y/^1^,1, 



(13) 



Eqs. (12,13) give the weak localization correction for any 
network. We now turn to few special cases. 
Incoherent networks. If all wires are longer than the 
phase coherence length (/^jy 3> i^, V (fJ-i^)), the weak 
localization correction (12) involves the same length as 
the classical conductance (10). If we write T^lp, — 
adNclel LcS-, then 

AT.,,, --^ {M-,')^^ = . (14) 

The ratio ATq,/ /T^]^, is network independent in this 
case. This simple result strongly relies on the hypothesis 



of wires with equal sections [14]. 

(b) Z„ 1 h 




FIG. 3. Two examples of mesoscopic devices. The wavy 
lines indicate connection to external reservoirs. 
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A ring. The transport through a ring has been studied in 
[16]. Here we rather consider the device on figure 3. a to 
illustrate the nonlocality of the weak localization. The 
classical conductance (10) is T|f] = ^xTlf^' inde- 
pendent of the presence of the arm and the ring. Then 
their weights are and this part of the network does not 
contribute to the weak localization correction (12). How- 
ever, since the cooperon Pc{x,x) is nonlocal, it feels the 
presence of the loop, even for x in the wires (3 1) 
and (1 ^4). Note that the naive uniform integration of 
the cooperon over the network (DRPM) strongly over- 
estimates the amplitude of the AAS oscillations (figure 
4) [14]. The decrease of the weak localization at high 
field (inset) is due to the contribution of the flux to the 
effective phase coherence length L^{(f)) [17]. 




FIG. 4. Dashed line : Act/ctq for a uniform integration of 
Pc- Continuous line : AT34/T'34 given by (12). The curves 
have been shifted so that they coincide at </> = 0. Inset : Same 
curves (without shift) for a higher window of flux (f>. The pa- 
rameters are la = h = ^ A*ni, Ic = 0.05 nm and Id = 5 /jm. 
W = 0.19 urn, L^{(t) = 0) = 1.7 urn. 

Multiterminal geometry. The origin of the negative weak 
localization ATa'0' lies in the negative weights, however 
for a multiterminal geometry some weights can be po.s- 
itive, like the weight(s) 9T^l^,/9/^^' of the wire(s) con- 
nected to other terminal(s) {e.g. fj! on figure 1). As a 
first example, we consider a wire on which is plugged 
one long arm of length 1^ connected to a third reservoir 
(figure 3.b for A^^ = 1). We focus ourselves on the fully 
coherent limit = oo. The classical (Drude) conduc- 
tance of this 3-terminal network is : T^l = ; ^^'^j^^l'.^i i ■ 
Then dT^\/dlc > 0. The wire [(2 ^ 1) + (1 ^ 3)] giveVa 
negative contribution to the weak localization correction 
whereas the arm (1 4) gives a positive one. Introduc- 
i'^g C//b//c = + + we find 

in the limit la//b ^ Ic- We now consider the case of Na 



long arms plugged in the middle of the wire {la = It) like 
on figure 3.b , to maximize their effect [14]. We obtain : 

AT23^i(-l + ^) , (15) 

a result valid for la <^ Ic ^ L^. We can now obtain a 
positive weak localization correction for Na > 4. This ef- 
fect is purely geometrical. Note that in the limit 1^ ^ 
the positive contribution vanishes. 

Conclusion. We have provided a general theory for the 
quantum transport of networks of diffusive wires con- 
nected to reservoirs. We obtained the classical conduc- 
tances and the weak localization corrections. We em- 
phasized the importance of the weights to give to the 
wires when integrating the cooperon over the network. 
This can lead to new geometrical effects like a positive 
correction (15), the physical reason being that coherent 
backscattering in a multiterminal geometry can enhance 
a transmission. 
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